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Abstract 

We studied a vertical "quantum dot molecule" , where one of the dots is occupied 
with electrons and the other with holes. We find that different phases occur in the 
ground state, depending on the carrier density and the interdot distance. When 
the system is dominated by shell structure, orbital degeneracies can be removed 
either by Hund's rule, or by Jahn- Teller deformation. Both mechanisms can lead to 
a maximum of the addition energy at mid-shell. At low densities and large interdot 
distances, bound electron-hole pairs are formed. 
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Semiconductor quantum dots share many qualities with natural atoms, and 
are thus often referred to as "artificial atoms" (for reviews see [1,2]). In single 
quantum dots, shell structure could be observed [3]: large differences in the 
addition energies occurred at particular electron numbers (N = 2, 6, 12) cor- 
responding to closed shells, forming the "quantum dot noble gases" . Just like 
in real atoms, for half-filled shells spin alignment due to Hund's rule leads to 
an increased stability of the system [3,4]. Quantum dots are uniquely suited 
to explore regimes of many-electron systems that are not accessible in conven- 
tional atomic physics, and much fundamental new insight on correlated states 
in these low-dimensional many-body systems was gained in the past years. 

Quantum dot artificial atoms can be coupled to form "artificial molecules". 
Such systems have been studied much both experimentally [5] and theoreti- 
cally [6], and were even discussed in connection with potential quantum com- 
puting devices [7] . Most work focused on vertically or laterally coupled double 
quantum dots with electrons. More recently, Anisimovas and Peeters [8,9] 
examined the correlated few-particle states in a vertical bipolar quantum dot 
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molecule, consisting of two vertically coupled quantum dots, one populated by 
electrons, the other by holes. They applied the exact diagonalization method 
in the limit of a strong, spin-polarizing magnetic field [8] and furthermore stud- 
ied the dynamic response within a hydrodynamic model [9]. The experimental 
realization of such electron-hole double quantum dots should for example be 
possible in bilayer-bipolar heterostructures with separated electron- and hole 
layers in equilibrium, as in biased GaAs/AlGaAs or InAs/GaSb [10]. These 
bipolar systems have attracted much attention recently, last but not least 
due to the exciting possibility of making a Bose-Einstein condensate with in- 
direct excitons [11]. The phases of a symmetric electron-hole bilayer system 
were recently also investigated by Liu et al. [12], and by Paolo et al. [13] in 
the fixed-node diffusion Monte Carlo approach [14] with periodic boundary 
conditions. 

In this Communication, we demonstrate that small bipolar-bilayer systems 
consisting of only a few electrons and holes show a very rich phase diagram 
which in fact appears to be different from the "bulk" [13]: depending on the 
electron- and hole-densities and the interdot distance, the ground state is 
either dominated by shell structure, leading to a competition between Jahn- 
Teller deformation [15] and spin alignment in the case of orbital degeneracies, 
or by the formation of localized, bound electron-hole pairs ("localized exci- 
tons"). 

As schematically sketched in the inset to Fig. 1, a quantum dot with electrons 
in one (two-dimensional) layer is separated by a distance zq from a quantum 
dot confining holes in the other 2D layer. The charge carriers in both layers 
are confined by a harmonic potential, V = m* eh ujQr 2 /2, with r 2 = x 2 + y 2 (a 
model which is frequently applied to describe single quantum dots [16,17]). 
The two layers are electrostatically coupled with no interdot tunneling, i.e., 
recombination is prohibited. 

To obtain the ground state energies and densities, we apply the density func- 
tional method in the local spin density approximation [18]. For the bipo- 
lar bilayer system, we have to solve four coupled equations for the electron- 
and hole spin densities, n|(r), n|(r) and n^(r),n^(r). For electrons and holes, 

respectively, the effective Kohn-Sham potential v^ a (r) consists of the ex- 
ternal harmonic confinement, the electron-electron (hole- hole) repulsion, the 
electron-hole (hole-electron) attraction and the exchange-correlation potential 
v a xc which is approximated in the local spin-density approximation. We use the 
von Barth and Hedin [19] formulation of the local electron (hole) exchange- 
correlation energy e xc {n e{h \ C [h) ) = e xc {n e{h \ Q)+f{C {h) )[e xc {n< h \ l)-e xc (n< h \ 0)], 
where is the electron (hole) spin polarization and f(( e ^) = ((1 + 

£e(ft))3/2 + (i _ £e(ft))3/2 _ 2 )/(2 3 / 2 - 2) is the polarization dependence. The 
function f(( e ^) interpolates between the paramagnetic (C e ^ = 0) and ferro- 
magnetic (C 6< ^ = 1) cases given by Tanatar and Ceperley [20]. The electron- 
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hole correlation is neglected. This is a good approximation for large interdot 
separation (the correlation energy is then diminished) and for strong exter- 
nal confinement (the correlation energy is then small compared to the kinetic 
and exchange energies). Even in the equilibrium density of the electron hole 
plasma the correlation energy is much smaller than the exchange energy [21] 
(the equilibrium density minimizes the total energy of the two-dimensional 
plasma in an infinite system). 

For example, the Kohn-Sham single particle wave functions for the electrons in 
one of the two dots, ^>f CT (r), satisfy the Kohn-Sham equations with the effective 
potential for the electrons, given by 



^(r) = im^ V + < c K(r),C e (r)) + J yf^dr' 

dr'. (1) 



n h (r') 



^\v-v'\ 2 + zl 
The electron density n e (r) in Eq. 1 is given by 

n e (r)=£l^(r)| 2 , (2) 
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where a corresponds to the two spins a = (|, t) an d i refers to the lowest 
occupied states. A corresponding relation holds for the hole density n h (r) in 
the other dot. C e ( r ) is the spin polarization for electrons, C e ( r ) — ( n f( r ) — 
nj(r))/n e (r). The Kohn-Sham equations for electrons and holes were solved 
self-consistently using a plane wave technique applied earlier to single dots [22,4] 
No symmetry restrictions for the charge or spin densities were applied. 

For simplicity, let us assume equal effective masses of electrons and holes, 
m* = m* h = m*, and also restrict the numbers of electrons and holes confined 
in each of the wells to be equal, N e = Nh = N. (Throughout the paper, we use 
effective atomic units with the energy measured in effective Hartree, Ha* = 
m*e 4 /ft 2 (47re e) 2 and the distance in effective Bohr radii, a* B = h 2 (4ire e)/m*e 2 . 
With the effective mass m* and dielectric constant e, the results then scale 
to the value of the semiconductor material in question.) The vertical double 
dot is then described by three parameters: the distance between the dots zq, 
the strength of the confining potential u (taken to be equal for both dots), 
and the number of electrons and holes. Due to the attractive interactions 
between the electrons and holes and since we assumed that N e = N h = N, 
the obtained ground state electron density necessarily equals that of the holes, 
n e (r) = n h (r); in the z-direction the densities are of course separated by zq. 
The spin density of the electrons and holes are also equal, but the relative 
direction of the electron and hole spins is not determined since there is no 
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direct spin-spin interaction. 

In the phase space determined by z and u three limiting cases are of partic- 
ular interest: 

(i) If both Zo and ujq are zero, the limit of the two-dimensional electron-hole 
plasma is reached. The average density for large TV will approach the equilib- 
rium density of the plasma. In the case of small droplets, for non-closed shells 
the orbital degeneracies are removed by internal Jahn- Teller deformation of 
the density, resulting in intriguing geometries of the ground state densities. 
This regime was called earlier the "ultimate jellium" limit, and has also been 
studied in connection with the physics of simple metal clusters [23] both in 
the two-dimensional [22] and in the three-dimensional case [24,25]. 

(ii) When ooq is still zero (or very small), but z$ becomes larger, the average 
particle density in each dot decreases and eventually, the particles begin to 
localize in a lattice. In comparison to the formation of "Wigner-molecules" in 
quantum dots with electrons [26] , here the localization is clearly seen also in the 
local density approximation due to the the attractive interaction between the 
electrons and the holes: they localize on top of each other forming a lattice of 
bound electron-hole pairs ("localized excitons"). Since the electrons and holes 
are confined in different layers, these electron-hole pairs have electric dipole 
moments and consequently repel each other. In each dot the localized electrons 
(or holes) show antiferromagnetic order (i.e., the so-called "spin density wave" 
[27,28,29]) in the ground state. The formation of bound electron-hole pairs is 
somewhat similar to the states discussed by Anisimovas and Peeters [8] for 
strong interdot coupling in a polarizing magnetic field. 

(iii) When ujq increases the strong confinement eventually hinders the internal 
deformation and the electron and hole densities are forced to be azimuthally 
symmetric. The degeneracy in an open-shell dot can then be reduced by mag- 
netization according to Hund's first rule [3,4]. 

Figure 1 shows the phase diagram of the N = 4 system. This case is particu- 
larly interesting since for a single dot, the second shell is half-filled having two 
p-electrons and spin S — 1 [4]. In the limit of the "ultimate jellium" (z = 0, 
ujq — 0), on the other hand, the system is strongly deformed and non-magnetic. 
Figure 1 shows the regions of different magnetic and geometrical structures 
in the zo-ujQ-plane. It is interesting to note that for small interdot distances 
the deformed system is the ground state up to a large value of u (for z = 
the circular S — 1 state becomes the ground state when ftw > 1.3 Ha*). For 
small values of u we see a transition from the deformed state to the state 
with localized excitons when z Q increases. 

A single dot with N = 6 has a filled electronic shell and naturally S = 0. 
Nevertheless, for a very weak confinement the ground state shows a spin- 
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Fig. 1. Phase diagram of an electron-hole double dot with 4 electrons and 4 holes, ujq 
is the strength of the external confinement and zq the interdot distance (in atomic 
units). The dashed lines show the estimated phase boundaries separating different 
spin-structures. The characteristic particle density is shown for each phase. The 
squares (spin S = 1) and triangles (S = 0) indicate the calculated points which 
are closest to the estimated phase boundaries. Inset to the right: Schematic picture 
of the bipolar quantum dot molecule with electrons (e) and holes (h) harmonically 
confined in two layers, which are separated by a distance zq. 



density wave [4]. A similar behavior is obtained for the electron-hole double 
dot at large values of the interdot distance zq. In the limit of small external 
confinement, the ground state appears to be a hexagon of localized excitons. 
Classically, a more favorable geometry would be a pentagon with one exciton 
at the center. However, in this case the hexagon is favored due to the formation 
of a spin-density wave, i.e., antiferromagnetic order of the spins in each dot (a 
centered pentagon would lead to "frustration"). The ground state of the N = 6 
double dot in the limit of the "ultimate jellium" (u;o = 0, z = 0) does not 
have azimuthal symmetry, but is triangular, which seems to be the dominating 
geometry for closed-shell electron- hole clusters in two dimensions [22,30]. The 
phase diagram (not printed here) shows a transition from a nonmagnetic trian- 
gular density to circular density when uq increases, and to localized electrons 
with a spin density wave when zq increases (for small u ). 
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Fig. 2. Minimum of the local r s as a function of the number N of electrons and 
holes for different values of the interdot distance zq. The average carrier density of 
a single dot was kept approximately independent of N (see text). 
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Fig. 3. The addition energy spectrum for different values of zq. The average carrier 
density of a single dot was kept approximately independent of N (see text). 
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In the above phase diagram, the number of electrons and holes N was fixed. 
When N is varied, however, experimental data on single quantum dots sug- 
gest [2] to keep the average density in each dot constant. This can be achieved 
by the relation u%(N) = e 2 /(47ree m*rf y/N) (see Ref. [4]). The parameter 
r s = Xj^fiva measures the in-layer coupling strength and is defined in terms 
of the in-layer areal density n e ^ h \r) for the electrons and holes, respectively. 
Here, we use the r s -value corresponding to the equilibrium density of the 
two-dimensional electron gas (r s = 1.51 a* B ). Neglecting the electron-hole cor- 
relation, this density is also the equilibrium density of the electron-hole plasma 
without external confinement. 

In the double dot system a single dot corresponds to the limit of infinite zq. 
When the interdot distance decreases, the electron-hole attraction increases 
the effective confinement and the density increases slightly. However, since the 
revalue has been chosen to be at the equilibrium density of the electron-hole 
plasma, r s = 1.5 a* B , this effect is small. Figure 2 shows the minimum local 
value of r s (determined from the maximum density) obtained in the dots for 
different values of N and z , using the above relation for u (N). For z = 3 a* B 
the minimum r s is already close to 1.5 a* B for all values of N. The variation as 
a function of N is caused by the shell structure: The smallest value is obtained 
for the closed shell system, N = 2, and the maximum value for the half-filled 
shell N = 4. The lowest curve (zo = 0) corresponds to the situation where 
the deformation wins against Hund's rule, leading to a different dependence 
of the minimum r s on N. 

In the case of single quantum dots the addition energy spectrum gives informa- 
tion about the shell structure. Maxima in the addition energy differences were 
observed for closed shells, and at mid-shell spin alignment due to Hund's rule 
leads to enhanced stability [3]. For the bipolar bilayer system studied here, 
we find that intriguingly, the addition energy spectra also reveal the internal 
Jahn- Teller deformation of the system. The difference in the electrochemical 
potentials of a bipolar double dot confining (N + 1) and N electrons and holes 
is given by A(N) = E(N+ 1) -2E(N) +E(N- 1), i.e., the second differences 
of the corresponding total ground state energies E(N). 

Figure 3 shows the addition energy differences A(JV) for z < 3 a* B . The addi- 
tion energy maxima for the closed shells at iV = 2 and N = 6 can be clearly 
seen. An additional maximum occurs at mid-shell, N — 4. It is interesting 
to note that also the case z = shows a similar spectrum. In this case the 
maximum at iV = 4 is not caused by Hund's rule, but by the Jahn- Teller 
deformation, which is strongest at half-filled shells. 

In conclusion, we have studied symmetric vertical electron-hole double dots 
using the geometrically unrestricted spin-density functional method. In the 
limit of small distances between the dots, the effective confinement is gov- 
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erned by the electron-hole electrostatic attraction and the properties of the 
system are determined by the Jahn- Teller deformation of the particle densi- 
ties. In the limit of weak external confinement, at increasing interdot distance 
electrons and holes form separate "excitons" which localize forming a lattice 
with antiferromagnetic spin-coupling in each dot. When the external confine- 
ment becomes large, Hund's rule determines the spin in each dot. The addition 
energy spectrum shows a maximum at half-shell. At large interdot distances 
this maximum is caused by Hund's rule while at small distances the maximum 
is a result of the Jahn- Teller deformation of the particle density. 

This work has been supported by the Academy of Finland under the Finnish 
Centre of Excellence Programme 2000-2005 (Project No. 44875, Nuclear and 
Condensed Matter Programme at JYFL), NORDITA, the Swedish Research 
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